Aging in Dense Colloids as Diffusion in the Logarithm of Time 
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The far-from-equilibrium dynamics of glassy systems share important phenomenological traits. 
A transition is generally observed from a time-homogeneous dynamical regime to an aging regime 
where physical changes occur intermittently and, on average, at a decreasing rate. It has been 
suggested that a global change of the independent time variable to its logarithm may render the 
aging dynamics homogeneous: for colloids, this entails diffusion but on a logarithmic time scale. Our 
novel analysis of experimental colloid data confirms that the mean square displacement grows linearly 
in time at low densities and shows that it grows linearly in the logarithm of time at high densities. 
Correspondingly, pairs of particles initially in close contact survive as pairs with a probability which 
decays exponentially in either time or its logarithm. The form of the Probability Density Function 
of the displacements shows that long-ranged spatial correlations are very long-lived in dense colloids. 
A phenomenological stochastic model is then introduced which relies on the growth and collapse 
of strongly correlated clusters ("dynamic heterogeneity"), and which reproduces the full spectrum 
of observed colloidal behaviors depending on the form assumed for the probability that a cluster 
collapses during a Monte Carlo update. In the limit where large clusters dominate, the collapse 
rate is oc 1/t, implying a homogeneous, log-Poissonian process that qualitatively reproduces the 
experimental results for dense colloids. Finally an analytical toy-model is discussed to elucidate 
the strong dependence of the simulation results on the integrability (or lack thereof) of the cluster 
collapse probability function. 

PACS numbers: 



I. INTRODUCTION 

Aging in amorphous materials has attracted 
widespread experimental, simulational and theoreti- 
cal interest for more than thirty years As a 
spontaneous off-equilibrium relaxation process, aging 
entails a decrease of the free energy and, correspond- 
ingly, a slow change of thermodynamic averages. E. g., 
in numerical studies of models for disordered mag- 
nets, the thermal energy decreases intermittently and, 
on average, at a decelerating rate during the aging 
process [6|-l9|]- Like these low-temperature materials, 
colloidal suspensions at high density (i. e. high volume 
fraction) exhibit intermittent dynamics and a gradual 
slowing down, here, in the rate at which particles move 
during light scattering [l(| [ll| and particle tracking 
experiments [L^ - [l5| . While the phenomenology of 
colloidal aging is broadly similar to that of thermally 
activated aging, spatially averaged quantities as energy 
and particle density hardly change in colloids. Further- 
more, no external field is required to elicit a measurable 
response. Finally, the motion of colloidal particles is 
time-homogeneous for sufficiently low density, but for no 
accessible time scales and for no value of the density does 
it look stationary or equilibrium-like. Generalizations of 
the Fluctuation Dissipation Theorem and the concept 
of effective temperature [lg, 1 17] , developed to describe 
thermally controlled aging in disordered magnets, are 
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not applicable to colloids, since the dichotomy between 
equilibrium-like fluctuations and off-equilibrium dynam- 
ics is absent. Reconciling differences and similarities 
between aging in colloids and other glassy systems thus 
requires a novel and broader approach . 

In this work, the colloidal particle tracking experi- 
ments of Courtland et al. [13] are re-analyzed, lead- 
ing to the conclusion that particle motion is diffusive 
in either time or logarithmic time, at low or high den- 
sities, respectively. Dynamics properties, such as the 
Mean Squared Displacement (MSD) and the persistence 
of tracer particles, appear homogeneous when viewed as 
a function of an appropriately re-scaled time variable. 
As the MSD, for instance, is a sum of independent dis- 
placements, the correct choice of time variable, hence, 
trivializes colloidal dynamics [30]. To explain why that 
is, a novel stochastic model is introduced with a sin- 
gle free parameter controlling the transition from diffu- 
sive to log-diffusive behavior. In the log-diffusive regime, 
the dynamics is driven by increasingly rare events, called 
"quakes", whose temporal distribution is homogeneous in 
logarithmic time. This links the dynamical behavior of 
dense colloids to the notion of "record dynamics" [HI, [lj| , 
a relaxation paradigm already applied to a number of 
complex systems d, |{§ l2CH2lj. 

Our experimental data analysis considers the MSD and 
the persistence, i. e. the probability that a pair of parti- 
cles initially close to each other remain close for at least 
time t, under a number of different conditions. The Prob- 
ability Density Function (PDF) of the particle displace- 
ment is likewise measured to illustrate intermittency and 
the long range of the spatial correlations. The model 
encodes the spatial heterogeneity in particle mobility, 
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Figure 1: (Color online) Left: Mean Square Displacement (in p 2 ) for colloids at low density, volume fractions p ~ 
0.46, 0.52, 0.53, 0.56, top to bottom, all plotted versus time, in seconds. The top-most data set has, for graphical reasons, 
been scaled down by a factor of five. Middle: same quantity for dense colloids, where for all data p ~ 0.62. Waiting times, 
in seconds: t w ~ 1, 10, 20, and 150 for pluses, squares, circles, and stars, respectively. Right: At high density, (p w 0.62) the 
persistence decays as a power law dependent on a separation-threshold 9. The inset shows the exponential decay at low density, 
p ~ 0.46, for the same thresholds. 



which is ubiquitous in glassy systems [TU, [TH, [U l23j . 
in terms of the formation and evolution of clusters of 
strongly correlated spatially contiguous particles. These 
clusters signal emerging "dynamic heterogeneities" (DH), 
i. e. a broad distribution of time-scales for the sur- 
vival and break-up of correlation patterns in different 
regions [12j. The stochastic model abstracts the under- 
lying dynamics solely into a cluster survival probability: 
complex behavior replaces simple diffusive behavior once 
clusters attain non-zero probability to survive indefinitely 
under the update rule. When the probability per up- 
date to break up falls off exponentially with cluster size, 
events in the model become well-separated in time, with 
a log-Poissonian statistic. This reproduces the qualita- 
tive features of the experimental data: the dynamics is 
intermittent and homogeneous in logarithmic time. 



II. EXPERIMENTAL DATA ANALYSIS 

Obtained using particle tracking techniques, the col- 
loidal data of Courtland et al. [13| comprise the 3d- 
trajectories of a set of tagged colloidal particles. The 
particle radius is 1.18 microns, and the position coordi- 
nates are given in microns, a length unit nearly equal to 
the particle radius. Colloids were initially centrifuged to 
obtain the desired density, and then briefly stirred. The 
end of the stirring phase is considered as the origin of 
the time axis, i.e., t — 0. Without time-translational 
invariance, this choice impacts the description of dense 
colloids. 

The MSD, relative to any overall drift in the system, 
is calculated from the trajectories as a function of time. 
Fig.[TJ leftmost panel, shows the well-known diffusive be- 
havior of low density colloids. The diffusion coefficient 
decreases monotonically with increasing volume fraction. 
MSD pertaining to high-density colloids are plotted in 
the middle panel on a log horizontal scale versus the scal- 



ing variable t/t wi where t w ,t w < t is the age at which 
data collection commences. Although a more precise 
knowledge of t w would facilitate a better data collapse, 
it is apparent that the particles in a dense colloid clearly 
diffuse in logarithmic time. 

The rightmost panel of Fig. Q] displays data for the 
persistence defined above. To ensure statistical indepen- 
dence, we partition the system in 20 x 20 x 5 identical 
sub- volumes. At time t wi a particle centrally located in 
each sub-volume and its closest neighbor are tagged as 
a pair. If at time t > t w the distance between the ele- 
ments of a pair increases beyond a given threshold value 
9 and subsequently remains above 9, the pair is deemed 
to have broken up at t. In the plot, three 9- values are 
considered. Measured in units of particle radius, these 
9 are large compared to the typical range of the relative 
motion. The persistence decays as a power of t for dense 
colloids, i. e. exponentially in log(i), while at low density 
it relaxes exponentially, as shown in the insert. 

Finally, Fig. [2] shows the empirical PDF (also known 
as the self-part of the Van Hove distribution function) 
of the absolute value of the particle displacement (in 
any direction) in the liquid (left panel) and glassy (right 
panel) regimes of the colloids. The vertical scale is log- 
arithmic. The three data sets displayed correspond to 
different values of the time interval At over which the 
displacement is observed. The inserts include more val- 
ues of Ai, and merely confirm the diffusive, respectively 
log-diffusive, character of the dynamics previously dis- 
cussed. Note that the MSD does not cross the origin. 
The particles already have acquired some dispersion af- 
ter the first sweep of the confocal microscope, probably 
due to a residual stirring motion. 

In similar analyses performed on a number of differ- 
ent systems [24l . [25j . a central Gaussian part is gener- 
ally flanked by an exponential 'non-Fickian' tail. The 
left-most part of the PDF should, by analogy, have a 
parabolic shape, and some curvature is indeed visible in 
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Figure 2: (Color online) The left panel pertains to the liquid, and the right one to the glassy regime of a colloid. Both panels 
display on a log-scale the PDF for the absolute value |Ar| of the particle displacement (in any direction). Three data sets are 
shown, corresponding to different durations (in seconds) of the time interval At over which the displacement is measured. In 
the inserts, the Mean Square Displacement (MSD) over a larger set of time intervals is plotted versus the interval length, using 
a linear (left panel) and logarithmic (right panel) abscissa. 



the left panel at small values of the abscissa. However, 
the exponential 'tail' extends right up to |Ar| = in 
the glassy regime. Thus, long-ranged spatial correlations 
survive over long times in dense colloids. This very fea- 
ture is incorporated as a central hypothesis in the model 
introduced next. 



III. CLUSTER MODEL 

The cluster model discussed below assumes that strong 
kinematic constraints bind the particles together inside 
clusters that form a spatially heterogeneous patchwork 
of particle properties. The growth and destruction of 
these clusters are the only mechanisms allowed for the 
net displacements of particles. The update dynamics is a 
Markov chain where a randomly selected cluster of (inte- 
ger) size h either survives intact or is destroyed with prob- 
ability P{h). By assumption, < P(h) < 1 decreases 
with h, i. e. larger clusters are more stable than smaller 
ones. The particles released by a collapsed cluster join 
neighboring clusters, and are activated to move in real 
space by a unit step in a random direction. Model vari- 
ants are characterized by different forms of P{h). Merely 
the integrability of P(h) in h discriminates between the 
behaviors similar to those of low and high density col- 
loids, and without loss of generality we can limit ourselves 
to 

P a (h) = — a = l,2..., (1) 

an a-family of models with P a=oa (h) = cxp(— h). 

In our Monte Carlo simulations the clusters are ar- 
ranged in a periodic lattice of length L. Initially, each lat- 
tice site holds uniformly a cluster of size hi — h. The av- 



erage number of particles bound in clusters Xh=i hi/L — 
h remains conserved. Note that these clusters do not 
per-se represent individual colloidal particles and their 
motion, since particles typically occupy space uniformly, 
with density fluctuations difficult to discern. Instead, 
ranges in cluster sizes represent DH in particle mobil- 
ities, which can shift dramatically even though actual 
particles hardly move. Here, heterogeneities arise dy- 
namically as larger clusters, once formed, are more long- 
lived than smaller ones, and the particles which belong 
to them correspondingly are less mobile. 

We update t — > t + 1 in parallel at all even and odd 
sites in alternating order. Each site i with hi > is 
updated by drawing a univariate random number r,. If 
ri < P a (hi), the cluster on that site breaks up into two 
randomly split fractions that are respectively added to 
the neighboring sites i ± 1, leaving i empty. To simulate 
the actual (activated) motion of particles, we add (mu- 
tually non-interacting) tracer particles which reside on 
lattice sites. Those walk randomly i — >• i± 1 iff there is a 
cluster on their site i that happens to break up. In par- 
ticular, a tracer particle can be stranded on a cluster-free 
site for a long time until a nearby cluster breaks up such 
that its scattered "debris" can reactivate that particle. 

Fig.[3]shows the MSD averaged over all tracer particles 
for a ~ 1 and a — oo. The former behaves diffusively (up 
to a cutoff scale t s t oc L 2 ). In the latter, no stationary 
state is approached and the behavior is diffusive in loga- 
rithmic time, as in the corresponding experimental result 
shown in Fig. [T] Fig. [3] also depicts persistence curves, 
i. e. probabilities that pairs of tracer particles initially 
located at the same grid point remain within a separa- 
tion of sites after t sweeps. As in the experiments, the 
(asymptotic) decay goes from being exponential in time 
for a — 1 (not shown [H[) to being exponential in log- 
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Figure 3: (Color online) Simulation results for different system sizes L for P 1 (h) ~ 1/h (left) and P°°(h) ~ exp(-ft) (center 
and right), to be compared with Fig. [T] For a — 1 (left) the MSD behaves diffusively, represented by the dotted line of unit 
slope, while for a — oo the MSD grows linearly in logarithmic time, with no sign of saturation. The persistence for a = oo 
(right) exhibits a power-law, independent of the threshold 9. 
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Figure 4: (Color online) Plot of the rate of cluster collapse vs. system age t for P (h) (left) and P°°(h) (right). 



arithmic time when the cluster collapse probability goes 
from hyperbolic to exponential at a — oo, independent of 
8. Unlike the experiments, changing 6 does not effect the 
power-law exponent but modifies the time for its onset. 



IV. RECORD DYNAMICS 

As the collapse of large clusters controls all aspects 
of the model in the aging regime, the statistics of such 
quake-events is key to the analysis, a situation identi- 
cal to aging in quenched glasses [7|. In turn, the fast 
collapse of small clusters provides a background of mo- 
bile particles diffusing within regions ("cages") bounded 
by large clusters. The divide between small and large 
clusters was empirically taken to be max(/i, <y(t)), where 
a(t) is the (growing) standard deviation of the cluster 
size distribution. The quake rate versus time, scaled by 
the system size L, is shown in Fig. 0] for our simula- 



tions. For P 1 (/i) the rate approaches a constant on a 
time scale independent of system size. For P°°(h), the 
stationary regime is in fact invisible even for the smallest 
systems, and the dynamics is thus purely non-stationary 
with a decelerating quake rate decaying as 1/t. Assum- 
ing that successive quakes are statistically independent 
events (demonstrated elsewhere [26]), they provide the 
true clock of the dynamics, and the particle MSD simply 
becomes proportional to the time integral of the quake 
rate: linear in time in the stationary (liquid) regime, and 
~ log(t) in the aging regime. Clearly then, the MSD in 
the aging regime between times t w and t scales as t/t w , a 
so-called full aging behavior which is consistent with the 
experimental findings. 

The model dynamics for the exponential col lapse p rob- 
ability can be linked to record dynamics [^ Il8l - l2ll[27l[28| . 
using the wide separation of time scales associated with 
the collapse of clusters of different sizes. Let us first 
note that the empirically estimated boundary between 
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small (i.e. irrelevant) and large (i.e. important) clusters, 
h(t) ks hit can in this case also be obtained by equat- 
ing to the time t the number of queries that a cluster of 
height h on average survives. The pace of the dynamics 
is always controlled by the smallest among the large clus- 
ters, henceforth the SL cluster. The latter is the first one 
to collapses with overwhelmingly high probability. Let 
the corresponding quake happen on a time scale . The 
new SL cluster can either be pre-existing or can be formed 
through the accretion process of the debris produced by 
the last quake, thus recapitulating all previous history, 
an accretion process which is now instantaneous on the 
time scale t^°\ Since the cluster size distribution is essen- 
tially continuous, the new SL cluster will only be slightly 
larger than its predecessor, a hallmark of "marginal sta- 
bility increase" Q. Hence, a record-sized random num- 
ber (corresponding to a record sized re-arrangement in 
the colloid) will suffice to destroy the new SL cluster. 
At time-step t the probability is 1 Jt that such event will 
happen. Statistical independence is assured by the wide 
separation of break-up events. Thus, inter-event times 
T («) = lni^-ln^™- 1 ) =log [t^/t (n ~^] become Poisso- 
nian distributed [27l . |29| . To ascertain the log-Poisson 
nature of the quake time distribution, we collect the 
quake times and show (elsewhere [2a|) that the se- 
ries (t^\t^ 2 \ . . . ,tW, . . .) of logarithmic waiting times 
t^ 1 ' is iid with the same exponential distribution. In 
contrast, when the cluster collapse probability becomes 
long-tailed, a < oo, break-up events overlap and long- 
lived clusters are deemphasized. 

V. MEAN-FIELD DESCRIPTION 

Finally, we provide a mean-field model of the cluster- 
ing dynamics that elucidates the preceding results. The 
fundamental quantity considered is the average number 
of clusters nk,t of binary size h = 2 k (k = 0, 1, 2, . . .) ex- 
tant at time t. For k > 0, those fc-clusters are immobile 
and break up with probability p^ ~ P(2 k ), but merely 

into two k — 1-clusters. Eq. ([T]) reduces to p% <~ (2 fe ) a . 
Integrability of pk is conferred via P(h)dh ~ pkd(2 k ) oc 
2 k pkdk, i. e. Ylk^ Pk < oo for a > 1. To write down a 
self-consistent equation for nk,t-> we stipulate that during 
each sweep the number of fc-clusters nk,t for all k > de- 
creases through break-ups by a factor 1—pk, yet, it also 
increases through the break-up of k + 1-clusters into two 
/c-clusters (with weight 2pk+i). Unit (k — 0)-clusters are 
indivisible but mobile. Larger /c-clusters further grow via 
accretion of those small mobile clusters, but (somewhat 
artificially) only in groups of 2 k with probability oc l/2 fe . 
During a sweep, all unit clusters attach somewhere, i. e. 
no,t+i is entirely unrelated to no,t- Hence, we obtain 

/ 1 s . , n k-l,t n 0,t ,„s 

n k ,t+i = fik,t (1 — Pk) + 2pk+in>k+i,t H ^F/y — ( 2 ) 

for k > 0, and no,j+i = 2p\n\ it - The last term describes 
growth through the accretion of all unit clusters no,t, 



equally shared between all clusters, N t = ^2^=0 n k',t- Its 
form ensures the conservation of ^2^=0^ n k',t(= hL), 
as used in our simulations above. 

Defining au,t = 2 k pi~nk,ti we rewrite Eq. ([2]) as 

— (ak,t+i — a*;,t) = afc+i,t — o-k.t + ak-i,t£-k-i%- (3) 
Pk 

Amazingly, our analysis finds that one of the factors, 
/Cfe = 1/ (2 k pk) or 7t = no, t /N t , is always sufficiently 
small asymptotically as to ignore the last term of Eq. . 
The remaining linear equation has stationary solutions 
(for t ~ t + 1 — > oo) entailing au+i ~ at ~ const, i. e. 
"fc ~ 1/ (2 k Pk)- Since the stationary cluster size distribu- 
tion must have a finite integral (after all, N t < L < oo), 
nk must decrease sufficiently, demanding pk 3> 2 _fc or 
a < 1. In this regime, 7t — > 7^o < 1 but ak-ilCk-i <C afc 
justifies dismissal of the last term in Eq. ([3]). Conversely, 
for 2 k pk C 1 (a > 1), the aging regime, numerics sug- 
gests that nk.t is highly localized at some increasing 
k = k(t): clusters grow and attain a typical size h(t). 
Unit-size clusters hardly ever occur, no,* = almost al- 
ways, and although /C^ grows exponentially, % as an- 
nihilates the last term. Rescaling r = p^t then turns 
Eq. ([5]) into a "wave" equation, [d T — dk] a(k, r) rj 0, 
with characteristic fc ~ t = pfei. Using pfe = and 
h ~ 2 fc , we extract the dominant growth- law for clus- 
ter sizes, ~ f 1 /", or for the event-rates (quakes) of 
cluster-size increases: 



The scaling reproduces for a — > 1 + and for a — » oo the 
simulation results in Fig. [3] (and, correspondingly, those 
for any a in-between |26|). While the equations do not 
describe the displacement of tracer particles explicitly, 
the discussion above implies that the quake rate sets the 
scale for spread in time for any observable, be it h(t) or 
MSD: logarithmic for a — >• oo and linear (with a cut-off) 
for a < 1. 



VI. CONCLUSIONS 

Experimental data are analyzed showing that colloidal 
motion is diffusive in either time or logarithmic time, de- 
pending on colloidal density. A model explaining these 
findings is presented, based on the survival probability of 
clusters, i. e. highly correlated groups of particles, rep- 
resenting dynamically heterogeneous regions of different 
mobility: a net displacement of a particle in the model 
requires the collapse of the cluster to which the particle 
belongs, i. e. a re-arrangement of (possibly large and in- 
creasing) correlation patterns. The survival probability 
of a cluster increases with its size. Once large clusters ob- 
tain the capacity to survive indefinitely, P(h) Cl//i for 
h oo, non-stationary, aging behavior arises, character- 
ized by intermittency and memory effects. Remarkably, 
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the t/t w scaling of the experimental data (full aging) is 
only achieved in the model in the limit of an exponential 
P(h). corresponding to the extreme value a = oo of its 
parameter. In that sense, tuning a parameter such as 
a, say, by changing the volume fraction p in a colloid, 
does not describe a sharp phase transition. But, presum- 
ably, the transition between both extremes occurs over a 
nearly unobservable interval of p. 
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